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I .  INTRODUCTION 


•■V 


The  US  Army  Missile  Command  (MICOM)  is  engaged  in  a  research  program  to 
develop  an  advanced  guidance  and  control  system  for  future  (1990* s)  Army 
modular  missiles.  Principal  investigators  have  defined  several  technical 
areas  in  which  government  and  contractor  personnel  are  contributing  to  the 
overall  program  objectives.  One  of  these  technical  areas  is  the  Disturbance 
Accommodating  Control  (DAC)  theory  and  the  design  of  discrete- time  disturbance- 
accommodating  controllers  for  discrete- time,  sampled-data  control  problems. 

The  DAC  method  of  design,  using  a  combination  of  waveform-mode  distrubance 
modeling  and  state- variable  control  techniques,  was  developed  by 
Dr.  C.  D.  Johnson  of  the  University  of  Alabama  in  Huntsville.  As  a  tool  for 
design  of  controllers,  the  DAC  approach  permits  three  primary  modes  of 
distrubance  accommodation:  (1)  cancellation  (absorption)  of  disturbance 
effects,  (2)  minimization  of  disturbance  effects,  or  (3)  constructive  utiliza¬ 
tion  of  the  disturbances  as  an  aid  in  acconplishing  the  primary  control  task. 
These  disturbance  accommodations  are  realized  in  addition  to  the  usual  control 
efforts  required  to  satisfy  system  performance  requirements  without 
disturbances. 

Hiis  report  describes  the  digital  computer  implementation  and  the  results 
obtained  from  the  design  of  a  discrete-time  disturbance-accommodating 
controller.  Itie  digital  conputer  inpleraentation  and  the  results  presented  in 
this  report  were  obtained  by  using  the  disturbance  cancellation  mode  of  the 
discrete-time  DAC  design  on  a  linear,  time- invariant  plant.  The  controllers 
obtained  via  this  design  method  are  termed  disturbance  absorption  controllers 
for  regulating  the  plant  state  to  zero. 

II.  DIGITAL  COMPUTER  IMPLEMENTATION 

Oliis  section  of  the  report  presents  the  details  of  the  discrete- time  DAC 
program  design  method.  The  first  section  presents  the  equations  of  the  total 
system  to  be  implemented.  However,  these  equations  may  be  considered  to  con¬ 
sist  of  three  groups :  (1)  continuous- time  disturbances,  (2)  continuous -time 

plant  (system),  and  (3)  discrete-time  DAC.  The  second  section  presents  the 
FORTRAN  program  implementation  of  this  DAC  design  method  for  execution  on  the 
PDP-11/34  digital  computer. 
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A.  THE  DAC  DESIGN  EXAMPLE 


The  discrete-time  DAC  design  program  example  to  be  considered  is  a 
regulator  problem  that  involves  a  first-order  controlled  plant  acted  upon  by  a 
constant  piecewise  disturbance.  The  complete  block  diagram  of  the  continuous¬ 
time  system  (plant)  with  the  continuous-time  DAC  design  is  as  shown  in 
Figure  1.  From  the  block  diagram,  the  following  equations  may  be  immediately 
written. 


For  the  piecewise  disturbance, 


,  , .  at 

w(t)  =  c  .e 

t 

For  the  first-order  plant  (system) , 
x(t)  -  a  x(t)  +  b  VinT)  +  /  w{t) 

x(t)  *  /*  x{t)  dt  +  x(o) 

o 

y{t)  =  c  x  (t) 

For  the  discrete-time  DAC  controller  design, 


y(nT)  =  y(t) 

TMPl  =  \edT  (a  -  d)  y  (nT)]  /  of  \(edT  -  eaT)\ 

Up  =  Ky(nT)/c 

Z(nT)  =  [|]  U(n  +  1)T] 

z(nT)  =  UnT)  -  TMPl 

Ud  =  {af/{b{a  -  J)]|  \(edT  -  eaT)/(eaT  -  l)J  z  (nT) 
U(nT )  =  U  +  Vd 


£,  [  (n  +  1)  T] 


U(nT) 


(1) 


(2) 

(3) 

(4) 


(5) 

(6) 

(7) 

(8) 
(9) 

(10) 

(ID 


(12) 


ORGANIZATION  OF  PLANT  AND  DISCRETE-TIME  OBSERVER 
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*DAC  Controller 


NOTE:  To  achieve  plant  deadbeat  response,  design  k  as; 
aT 

k  -  ae 

db  - 7T 

b(l-e  ) 

Figure  1.  Organization  and  block  diagram  of  problem. 
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where  K 


(13) 


=  ae  / 


aT/  b( 1 


E  =  —  *  a  unit  delayer  (14) 

Observe  that  there  are  three  cases  to  be  considered  in  the  use  of 
Eqs .  (6),  (10),  and  (12).  These  three  cases  are  as  follows:  (1)  where  d  -  0, 

(2)  where  d  =  a,  and  (3)  where  neither  case  (1)  or  (2)  applies.  These  cases 
are  considered  in  the  DAC  program  implementation  by  using  the  FORTRAN  IF 
statement. 

Note  that  the  disturbances  W{t)  have  been  experimentally  modeled  and 
found  to  obey  a  first-order  equation  of 


z  =  dz  +  o  (t)  (16) 

where  the  (real)  coefficient  d  is  assumed  to  be  known.  Now,  the  disturbance 
model  (given  by  Eq.  (1))  is 


where  a  is  any  arbitrary  (real)  scalar  constant.  This  implies  that 


z(t)  =  c.e* 

t 


z{t)  =  a c  .e 

v 


=  dz(t) 


Therefore,  a  =  d,  which  may  be  obtained  from  these  relationships  as  shown. 

By  setting  d  -  0,  the  case  of  constant  disturbances  is  obtained;  setting 
d  >  0  yields  the  case  of  exponentially  growing  disturbances;  and  setting  d  <  0 
yields  the  case  of  exponentially  decaying  disturbances.  Thus,  this  permits  a 
wide  range  of  realistic  disturbances  to  be  utilized  in  the  design  process  of 
high  performance  digital  controllers. 

Also,  the  above  equations  have  been  written  in  the  mathematical  notation 

i  ^ 

of  Reference  1  as  utilized  by  Dr.  Johnson.  Table  1  is  shown  as  an  aid  to  — 

understanding  the  notation  of  the  above  equations  and  their  usage  in  the 
FORTRAN  DAC  program  implementation;  this  implementation  is  described  in  the 
next  section. 


TABLE  1.  A  CROSS  INDEX  OF  VARIABLES 


MATHEMATICAL 

SYMBOL 

FORTRAN 

SYMBOL 

VALUE  USED 

IN 

CASE  1 

MEANING  AND/OR  USAGE 

A,  a 

A 

1.0 

Coefficient  of  x(t) 

a 

AWT 

0.0 

Coefficient  of  ea 

B,  b 

B 

1.0 

Coefficient  of  plant  equation 

C 

C 

1.0 

Coefficient  of  plant  output 

C. 

t 

CWT 

1.0 

Coefficient  of  Ce0, 

M 

D 

D 

Calculated 

Same  as  AWT 

dt 

DT 

1/64 

Integration  interval 

f 

F 

1.0 

Coefficient  of  disturbance 

— 

ICN 

1.0 

Case  number  parameter 

— 

ITP 

0.0 

Print  interval 

k 

K 

Calculated 

Equation  (13) 

— 

KUTTA 

Calculated 

Integration  control  parameter 

— 

NX 

1.0 

Number  of  variables  to  integrate 

T 

ST 

8  •  DT 

Sample  interval 

t 

T 

0.0 

Time 

— 

TMPl 

Calculated 

Equation  (6) 

~ 

TSTOP 

1.0 

Time  to  stop 

ud 

UD 

Calculated 

Equation  (10) 

U(nT)  =  £/(i) 

UNT 

Calculated 

Equation  (11) 

V 

p 

UP 

Calculated 

Equation  (7) 

Wit) 

WT 

Calculated 

Equation  ( 1) 

x  it) 

XDT 

Calculated 

Equation  (2) 

U(«  +  1)7] 

XINPT 

Calculated 

Equation  (12) 

C  inT) 

XINT 

Calculated 

Equation  (8) 

xit) 

XT 

Calculated 

Equation  (3) 

yinT) 

YNT 

Calculated 

Equation  (5) 

yit) 

YT 

Calculated 

Equation  (4) 

zinT) 

ZHNT 

Calculated 

Equation  (9) 

— 

V 

Calculated 

Array  of  graphical  plot  variables 

1/E 

— 

— 

Used  to  denote  a  unit  delayer 

1 

I _ £ _ 

EXP 

Natural  logarithm  as  in  Eq.  (1) 

T 
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B.  THE  PDP  11/34  PROGRAM (S) 

This  section  of  the  report  describes  the  digital  computer  implementation 
of  the  first-order  system  (plant)  with  a  constant  piecewise  disturbance  and 
the  associated  designed  discrete-time  DAC  controller.  These  programs  have 
been  programmed  in  the  FORTRAN  programming  language  for  execution  by  a 
PDP  11/34  digital  computer  and  its  operating  system  software.  A  listing  of 
the  main  program  is  as  shown  in  Figure  2,  which  contains  the  equations  given 
in  the  previous  section  of  the  report  for  the  constant  piecewise  disturbance 
model,  the  first-order  plant  (system)  model,  and  the  discrete- time  DAC 
controller  design  model.  Also  contained  in  this  listing  are  the  equations  and 
necessary  statements  for  the  initial  conditions,  values  of  variable  parameters, 
and  control  of  the  program  during  execution. 

The  required  plant  differential  equations  (one  in  this  case)  are 
integrated  by  the  fourth-order  Runga-Kutta  integration  scheme  as  contained  in 
its  listing  shown  in  Figure  3.  This  very  familiar  and  widely  used  integration 


scheme  is 

^n+1 

-  yn  +  Hko +  2ki 

+  2 k2  +  fc3) 

(20) 

where  k 

=  h  •  -Kv  0 

(21) 

•  h  ‘  f(xn  +  I  h' 

*n  +  I  ko) 

(22) 

*2 

=  h  •  f(xn  *  i  h. 

y  +  t  k*) 

°n  2  1/ 

(23) 

*3 

=  h  •  /(*„  +  h,  yn  +  k2) 

(24) 

Table  1  is  a  cross  index  of  the  variables  contained  in  the  block  diagram 
(Figure  1  and  the  mathematical  equations  written  from  the  block  diagram)  and 
the  FORTRAN  program  variable  names  contained  in  the  listings  of  the  designed 
discrete- time  DAC  digital  conputer  implementation.  Also,  Table  1  presents  the 
value  of  the  variables  as  used  in  the  execution  of  the  DAC  program  for  case  1. 

The  graphical  plots  of  the  dependent  variables  versus  the  independent 
variable  time  for  the  designed  DAC  controller  program  were  obtained  by  the 
execution  of  the  plot  program  as  shown  in  the  listing  of  Figure  4.  A  legend 
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C  ***  DAC  PROGRAM  EXAMPLE  NUMBER  1* 

C  ***  FIRST-ORDER  SYSTEM  WITH  A  CONSTANT  DISTURBANCE . 

C  ***  UNSTABLE  SYSYTEM 
REAL  K 

COMMON  KUTTA f  DTf  NX*  XT»  XDTf  YT 
DIMENSION  0(12) 

DATA  A/1 *0/ 

DATA  B/l * 0/ 

DATA  C/1.0/ 

DATA  F/1.0/ 

DATA  UNT/O.O/ 

NX  *  1 
T  =  0.0 
DT  *  1.0/64.0 
TSTOP  =1.0 
ITP  •  0 
ST  =  8.0*DT 
XINPT  =  0.0 
PRINT  1021 
PRINT  20 
READ(5f21 )  ICN 
WRITE<2)  ICN 
PRINT  23 
READ ( 5  f 1040 )  XT 
PRINT  80 

REAEK5  f  1040 )  CWT 
PRINT  81 

READ<5f 1040)  AWT 
D  *  AWT 
PRINT  1020 
9  CONTINUE 

DO  51  KUTTA  =  1f4 
YT  ■  C*XT 

GO  TO  <  60f50f  30f40)  KUTTA 
60  CONTINUE 

IF(MOD( ITP?8>  .NE.  0)  GO  TO  62 

0(1)  -  XT 

0(2)  -  YNT 

V ( 3 )  *  UNT 

V ( 4 )  *  UD 

0(5)  -  UP 

0(6)  *  XINT 

0(7)  =  XINPT 

0(8)  =  TMP1 

0(9)  *  ZHNT 

0(10)  -  XDT 

0(11)  *  WT 

0(12)  «  T 

WRITE (2 )  0 

K  *  A*EXP(A*ST)/(B*<1.0  -  EXP(A*ST>>> 

^NT  ■  YT 

IF( D  .EQ.  0.0)  TMP1  =  A*YNT/(C*F*( 1 .0-EXP( A*ST) ) ) 

IF  (D  .EQ.  A)  TMP1  =  - ( (EXP(D*ST)*YNT)/(C*F*ST*EXP< A*ST> > ) 

IF < D  .NE.  0.0  .AND.  D  .NE.  A)  TMP1  «  EXP ( D*ST ) * ( A-D > *YNT / ( C*F* ( EXP 
1(D*ST>  -  EXP < A#ST ) ) ) 

UP  -  ( K/C  >  *YNT 

Figure  2.  Listing  of  the  PAACP1  main  program. 
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0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 


0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 


XINT  *  XINPT 

ZHNT  -  XINT  -  TMP1 

IF ( D  .EQ.  0*0)  UD  =■  ~(F/B)*ZHNT 

IF  ( D  .EG.  A)  UD  =  ( A*F*ST*EXF‘ < A*ST ) *ZHNT )/( B* ( 1  ♦  0-EXP ( A*ST ))  ) 

IF ( D  .NE ♦  0*0  .AND*  D  ♦ NE .  A)  UD  «  ( A*F* ( EXP ( D*5T )  -  EXP ( A*ST > ) *ZH 
1NT>/(B*(A-D)*<EXP(A*ST>-1.0>) 

UNT  =  UP  +  UD 

D  IF ( T  *LT .  ST)  UNT  *  UP 

IF ( D  .EG *  0.0)  XINPT  *  EXP ( A*ST ) *TMP1  -  (B/F)*UNT 
IF ( D  .EG.  A)  XINPT  *  - ( EXP ( ( A+D )*ST > *YNT ) / ( C*F*ST*EXP ( A*ST ) )  - 
1  <B*EXP(D*ST)*(EXP< A*ST>  -  1 ♦ 0 ) *UNT ) / ( A*F*ST*EXP ( A*ST ) ) 

IF ( D  .NE.  0.0  .AND.  D  .NE.  A)  XINPT  *  < < A-D> *EXP ( ( A+D ) *ST )*YNT ) / < C 
1*F*( ( EXP  (  D*ST ) -EXP  <  A*ST ) ) ) )  +  ( B* ( A~D ) *EXP ( D*ST ) *(EXP ( A*ST ) -1 ♦ 0 ) *U 
2NT ) / <  A*F*EXP ( D*ST ) -EXP ( A*ST ) ) 

62  CONTINUE 

WT  a  CWT*EXP(AWT*T) 

XDT  *  A*XT  -F  B*UNT  +  F*WT 

IF ( HOD < I TP *  04 )  .NE.  0)  GO  TO  61 

PRINT  1010f  T t  XDT *  XT*  YNT 

PRINT  1010*  UP*  UD*  UNT *  TMP1 

PRINT  1010f  XINPT ,  XINT  *  ZHNT *  N*  WT 

PRINT  1021 

61  CONTINUE 
Y(l)  *  XT 
V(2)  *  YNT 
0(3)  *  UNT 
0(4)  «  UD 
0(5)  *  UP 
0(6)  *  XINT 
0(7)  *  XINPT 
0(8)  a  TMP1 
0(9)  »  ZHNT 
0(10)  *  XDT 
0(11)  =  WT 
0(12)  a  T 
WRITE ( 2 )  0 

30  T  ®  T  +  0.50*DT 

40  CONTINUE 

50  CALL  RUNGK 

51  CONTINUE 

ITP  a  I TP  +  1 
IF ( T  .LE.  TSTOP )  GO  TO  9 
PRINT  1030*  DT*  ST*  WT*  XT 
STOP 

20  FORMAT (/*T12*  'DAC  PROGRAM  *1*  CASE  ♦'*♦) 

21  FORMAT (12) 

23  FORMAT ( T1 2  * ' INPUT  XT  a  '*♦) 

80  F0RMAT(T12* 'FOR  EXPONENTIAL  DISTURBANCES )  *  '  */* 

1T12* 'INPUT  CWT  *  '*♦) 

81  FORMAT (T 12* 'INPUT  AWT  a  ',$) 

1010  FORMAT (5(4X*E12.5) ) 

1020  FORMAT (/*/*T12* 'DAC  PROGRAM  EXAMPLE  NUMBER  l.'*/*T12* 

1 'OUTPUT  FORMAT?'*/*/* 

28X* 'TIME' * 12X* 'XDT' * 13X* 'XT  a  YT ' * 9X * ' YNT ' * /* 

38X * 'UP  ' * 12X * ' UD  ' * 12X* 'UNT  ' * 12X* 'TMP1 ' */* 

48X  * 'XINPT' *1 IX* 'XINT' * 12X * ' ZHNT ' * 12X • • K  ' *  1 2X * ' WT ' * / ) 

Figure  2.  Listing  of  the  DAAC^l  main  nrogram  (Cont’d). 
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nv  \Pi&  '  "  «- 


0099  1021  FORMAT ( 1H  ) 

0100  1030  FORMAT (/»  T12r'CASE  PARAMETERS: ' r/rT12r 

1 ' INTEGRATION  SCHEME*  RUNGA-KUTTA  4TH  ORDER' r/rT12r 
2 'INTEGRATION  STEP  SIZE:  DT  *  ' r E12.5r /rT12r 
3 'SAMPLE  INTERVAL:  ST  =  ' rE12.5r/r T12 
4 'DISTURBANCE:  WT  =  ' rE12.5r/rT12r 
5 'EQUATION  FOR  UNT:  UNT  *  UP  t  UD'r/rT12r 
6 'STEADY  STATE  OUTPUT:  X(T)  =  ' fE12.5f /»/»/) 

0101  1040  F0RMAT<F10*4) 

0102  END 


PROGRAM  SECTIONS 


NAME 

SIZE 

ATTRIBUTES 

♦C0DE1 

003330 

876 

RWr I ? CONrLCL 

*P DATA 

000012 

5 

RWr DrCONfLCL 

•IDATA 

001050 

276 

RUrDrCONrLCL 

•VARS 

000204 

66 

RU* DfCONfLCL 

♦TEMPS 

000010 

4 

RWfDfCONfLCL 

.♦••«. 

000024 

10 

RWfDfOVRfGBL 

TOTAL  SPACE  ALLOCATED  ■  004652  1237 


Figure  2.  Listing  of  the  DAACP1  main  program  (Concluded). 
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0001 

SUBROUTINE  RUNGK 

0002 

COMMON  KUTTAr  DT,  NX,  X,  DX 

0003 

DIMENSION  X(l>»  DX ( 1 > »  XA(1>»  DXA(l) 

0004 

GO  TO  (10r30r50r70) fKUTTA 

0005 

10 

DO  20  I  *  1» NX 

0006 

XA( I )  =  X(I) 

0007 

DXA(I)  =  DT*DX(I) 

0008 

20 

X(I)  =  X(I)  +  0.5*DXA(I> 

0009 

RETURN 

0010 

30 

TDT  =  2.0#DT 

0011 

HDT  =  0.5*DT 

0012 

DO  40  I  =  1, NX 

0013 

DXA(I)  »  DXA(I)  +  TDT*DX(I) 

0014 

40 

X(I>  a  XA( I )  +  HDT*DX<I) 

0015 

RETURN 

0016 

50 

DO  60  I  =  1 f NX 

0017 

VDT  =  DT*DX(I) 

0018 

DXA ( I >  a  DXA(I)  +  2.0*VDT 

0019 

60 

X(I)  =  XA ( I )  +  VDT 

0020 

RETURN 

0021 

70 

DO  80  I  =  1 *NX 

0022 

80 

X(I>  =  XA ( I )  +  (DXA(I)  +  DT*DX( I ) >/6.0 

0023 

RETURN 

0024 

END 

PROGRAM  SECTIONS 


NAME 

SIZE 

ATTRIBUTES 

♦C0DE1 

000510 

164 

RW» I »CON»LCL 

♦PDATA 

000012 

5 

RWrDfCONrLCL 

♦VARS 

000026 

11 

RU?D»CON»LCL 

,♦♦♦♦» 

000020 

8 

RUrDrOVRrGBL 

TOTAL  SPACE  ALLOCATED  *  000570  188 


Figure  3.  Listing  of  Runga-Kutta  integration  routine. 
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V* 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

oooa 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 


C  ***  GENERALIZED  PLOTTING  PROGRAM  -  FOR  DACS. 

DIMENSION  V < 12 ) 

DIMENSION  PT ( 1025 ) 

DIMENSION  PXT < 1025 ) 

INTEGER*2  DRAW  9  LOOP 

REWIND  2 

READ  (2)  ICN 

PRINT  22 1  ICN 

PRINT  100 

PRINT  32 

PRINT  23 

READ(5f21 )  DRAW 

IF (DRAW  «EQ*  7  Y7 )  CALL  HDCOPY 

CALL  ST7611 

CALL  INITT ( 240 ) 

DO  20  J  =  1» 12 
REWIND  2 
PRINT  30 
READ  <5f 31 )  NOP 
READ  (2)  ICN 
CALL  NEWPAG 
PRINT  122»  NOPr  ICN 
DO  10  I  3  If  1025 
READ ( 2  r END-66 )  V 
PT ( I )  *  0(12) 

PXT ( I )  =  O(NOP) 

NP  s  I 

10  CONTINUE 

66  CONTINUE 
CALL  BINITT 
CALL  NPTS(NP) 

CALL  SYMBL(l) 

CALL  SIZES (0,50) 

CALL  CHECK (PT f PXT) 

CALL  DSPLAY(PTfPXT) 

CALL  MOOABS ( 1 00 1 50 ) 

CALL  ANMODE 
PRINT  23 

READ  (5f21)  DRAW 

IF (DRAW  *EQ*  7Y7>  CALL  HDCOPY 

CALL  NEWPAG 

PRINT  33 

READ(5f21)  LOOP 

IF (LOOP  *EQ*  7Y7)  GO  TO  20 

CALL  HT7611 

REWIND  2 

STOP 

20  CONTINUE 

21  FORMAT (A1 ) 

22  FORMAT ( T1 2 9  7 DAC  PROGRAM  *1*  CASE  #'»I2) 

122  FORMAT ( T15  9  7  PLOT  NO*  79l2r7*  DAC  PROGRAM  41 9  CASE  ♦/9l2f7*7) 

23  FORMAT ( 1 X  9  7  * 7  9  ♦ ) 

30  FORMAT </ 7/9 IX 9 'ENTER  *  OF  PLOT(S)  WANTED  (UP  TO  12)  'til 

31  FORMAT (13) 

32  FORMAT ( 1 X 9  7  AFTER  STAR  APPEARS  ON  THE  DISPLAY 7 9 /w IX 9 

1 'ENTER  A  ♦ Y ♦  AND  A  *CR*  TO  DRAW  A  HARDCOPY  OF  THE  DISPLAY*7) 

Figure  4.  Listing  of  the  graphical  plot  program. 
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0055 

0056 


0057 


33  FORMAT </,/,' 
100  FORMAT </»/»Tl 

*' - 


1 'PLOT 
2 'PLOT 
3' PLOT 
4 'PLOT 
5 'PLOT 
6 'PLOT 
7 'PLOT 
8 'PLOT 
9  'PLOT 
A 'PLOT 
B'PLOT 
C 'PLOT 
END 


NO. 

NO. 

NO. 

NO. 

NO. 

NO. 

NO. 

NO. 

NO. 

NO. 

NO. 

NO. 


1 t 
21 
31 
41 
51 
6! 
71 
81 
91 
101 
111 
121 


ANY  MORE  GRAPHS?  ',♦) 

2? 'LEGEND  TO  THE  PLOT(S) 1 ' »/,T12» 
- ' 

-  X(T>  VERSUS  TIME. ' >/»T12, 

-  Y(NT>  VERSUS  TIME. ' »/» T12r 

-  U(NT)  VERSUS  TIME. ' »/,T12* 

-  UD  VERSUS  TIME. ' »/»T12» 

-  UP  VERSUS  TIME. ' f/»T12» 

-  XINT  VERSUS  TIME. ' »/»T12» 

-  XINPT  VERSUS  TIME. ' f/»T12» 

-  TMP1  VERSUS  TIME. ' f/fT12» 

-  ZHNT  VERSUS  TIME. ' »/»T12» 

-  XDT  VERSUS  TIME . ' r/?T12» 

-  W(T)  VERSUS  TIME. ' »/»T12» 

-  TIME  VERSUS  TIME. '»/»/»/) 


PROGRAM  SECTIONS 


NAME  SIZE 


♦C0DE1 

001152 

309 

♦PDATA 

000024 

10 

♦I DATA 

001452 

405 

♦VARS 

020106 

4131 

ATTRIBUTES 

RWfIfCONfLCL 

RU»D,CON»LCL 

RUfDfCONfLCL 

RW»D»CONfLCL 


TOTAL  SPACE  ALLOCATED  =  022756  4855 


Figure  4.  Listing  of  the  graphical  plot  program  (Concluded) 


to  the  graphical  plots  capable  of  being  output  by  this  plot  program  from  the 
DAC  controller  program  is  shown  in  Table  2. 

C.  PROGRAM  EXECUTION 

This  section  of  the  report  presents  the  necessary  operating  system  software 
control  commands  for  the  PDP-11/34  digital  conputer  to  execute  the  designed  DAC 
controller  program  implementation.  Before  the  program — any  program,  for  that 
matter — can  be  executed  by  the  PDP-11/34  conputer,  a  task  for  the  program  must 
be  created.  Therefore,  included  in  this  section  will  be  the  necessary  control 
commands  to  create  the  task  of  the  discrete-time  DAC  controller  program 
described  in  the  previous  sections  of  this  report. 

The  following  assumptions  are  assumed  concerning  these  control  commands: 

1.  The  operating  system  has  been  booted  in  the  computer. 

2.  The  proper  user  identification  code  (UIC)  has  been  set. 

3.  The  program(s)  exists  on  the  disk  pack  in  source  code.  (The  program 
name  is  DDACP1.FTN. ) 

Then  the  control  commands  contained  in  Table  3  may  be  input  to  the  operating 
system  by  the  terminal  operator  onto  which  the  operating  system  is  logged  in 
to  create  the  task  for  the  DDACP1.TSK  program. 

To  execute  the  task  of  DDACP1.TSK  at  the  present  session,  or  any  future 
session,  the  operator  must  input  the  control  commands  shown  in  Table  4  onto 
the  input  terminal. 

A  similar  list  of  control  commands  must  also  be  input  for  the  graphical 
plot  program.  However,  this  source  and  task  program  is  given  the  file  name 
DACTKP .  It  is  further  assumed  that  the  DACTKP  task  is  logged  in  by  the  user 
onto  the  graphics  display  terminal  (device  TT1:) .  This  list  of  commands  is 
as  shown  in  Table  5. 

III.  RESULTS  OF  EXECUTION 

This  section  of  the  report  presents  the  results  of  the  DAC  design  example 
implementation  and  its  execution  by  the  PDP-11/34  digital  computer.  The  out¬ 
put  resulting  from  the  FORTRAN  READ  and  WRITE  statements  are  as  contained  in 
Figure  5.  These  results  are  also  contained  in  the  plots  of  the  variable  ver¬ 
sus  time  as  shown  in  Figures  6  through  17.  The  significance  of  the  symbol  0 
on  the  plots  (graphs)  is  that  the  symbol  appears  at  every  DT  seconds;  that  is, 
when  the  value  of  Kutta  •  1. 


TABLE  2.  LEGEND  TO  THE  GRAPHICAL  PLOTS 


LEGEND  TO  THE  PL0T(S)l 


PLOT  NO. 
PLOT  NO. 
PLOT  NO. 
PLOT  NO. 
PLOT  NO. 
PLOT  NO. 
PLOT  NO. 
PLOT  NO. 
PLOT  NO. 
PLOT  NO. 
PLOT  NO. 
PLOT  NO. 


1*  -  X(T) 
2*  -  V(NT) 
31  -  U(NT) 
41  -  UD 
51  -  UP 
6*  -  X1NT 
7*  -  XINPT 
8*  -  THP1 
91  -  ZHHT 
l#t  -  XDT 
11*  -  U(T) 
12*  -  TINE 


VERSUS  TIRE. 
VERSUS  TIRE. 
VERSUS  TIRE. 
VERSUS  TINE. 
VERSUS  TIRE. 
VERSUS  TIRE. 
VERSUS  TIRE. 
VERSUS  TINE. 
VERSUS  TIRE. 
VERSUS  TIRE. 
VERSUS  TIRE. 
VERSUS  TINE. 
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TABLE  3.  A  TYPICAL  LIST  OF  TASK  CREATION  COMMANDS 


F4P  DDACP1 . OBJ  -  DDACP1.FTN 
TKB  DDACP1.TSK  -  DDACPl.OBJ 
PIP  DDACP1.*;*/PU 
PIP  DDACP1 . OBJ; */DE 
PIP  DDACP1.*;*/LI 


TABLE  4.  A  TYPICAL  LIST  OF  TASK  EXECUTION  COMMANDS 

INS  DDACP1 
LUN  DDACPl 
REA  DDACPl  5  TL; 

REA  DDACPl  6  TI: 

LUN  DDACPl 
RUN  DDACPl 
PIP  FOR002 ,*;*/PU 
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TABLE  5.  GRAPHICAL  PLOT  PROGRAM  TASK  CREATION  COMMANDS 


>  F4P  DACTKP* OBJ • DftCTKP . FTN,NEUPAG . FTN 
^TKB 

TKB>  DACTKP .  TSK-DACTICP .  OB  J#  C277, 43AG2 .  OLB/LB 
TKB>/ 

ENTER  OPTIONS* 

TKB>ASG  -  TTlll 
TKB>ASG  •  TT1*\*\*6 
TKB>ASG  «  TT1*9 
TKB>// 

> 


BAG  PROGRAM  tl  r  CASF  tj. 

INPUT  XT  =  1.0 

FOR  EXPONENTIAL  MSTI IRPANCE < 8  )  l 
INPUT  CUT  ™  1.0 
INPUT  AUT  =  0.0 


DAC  PROGRAM  EXAMPLE  NUMBER  1. 

OUTPUT  format: 

TIME  XDT  XT  “  YT  YNT 


UP 

XINPT 

UD 

XINT 

UNT 

ZHNT 

TMF‘l 

K 

UT 

0 ♦ OOOOOE+OO 
-0 * 85104E+0.1 

0  *  75 104E+0 1 

-0. 14021E+02 
-0*  75104E+01 

0 ♦ OOOOOE+OO 

0* 10000E+01 
-0* 1A021E+02 

0 ♦ 75104E+O1 

0* 10000E+01 
-0  *  75104E+01 
-0 ♦ 85104E+01 

0 . 10000E+01 

0 • A2500E-01 
-0.85104E+01 

0 ♦ 75104E+0I 

-0* 14918E+02 
-0  ♦  75104E+01 

0 . OOOOOE+OO 

0* 10294E+00 
-0  ♦  1 A021E+02 

0 ♦ 75104E+01 

0  * 10000E+01 
-0 ♦ 75104E+0L 
-0 ♦ 85104E+0 1 

0. 10000E+01 

0* 12500E+00 

0. 72447F+01 
0*11153E+0:i 

0 ♦ 62799E+01 
-0 . 1 1 153E+01 
0*75104E+01 

-0  *  85150E+00 

0* A1314E+01 
0*1 1 153E+01 

-0  *  85150E+00 

0  *  A395IE+0 1 
-0  *  85104E+01 

0 . 10000E+01 

0* 1 8750E+00 
0.724A7E+01 

0*  U153E+01 

0 ♦ 668 1 7E  +  01 
-0*11153E+01 

0 ♦ 75 1 04E+01 

-0.44972E+00 

0 ♦ 61314E+01 

0* 11153E+01 

-0 ♦ 851 50E+00 

C. A3951E+01 
-0 ♦ 85104E+01 

0.10000E+01 

0 ♦ 25000E+00 

0. 189l7Ff00 
0.94837E+ 00 

0 ♦ 21853E+00 
-0 . 94837E+00 
0*11 153E+01 

-0 ♦ 22222E-01 
-0 ♦ 75925E+0O 
0.94837E+00 

-0  *  22222E-01 

0 ♦ 1A690E+00 
-0  *  85104E+01 

0. 10000E+01 

0 ♦ 31250E+00 

0, 19912E+00 

0  *  94837E+00 

0 » 2325 1E+00 
-0 ♦ 94837E+00 

0* 1 1153E+01 

-0*82408E~02 
-0 ♦ 75925E+00 

0 * 94837E+00 

-0*  22222E-01 

0 ♦ 1&690E+00 
-0. 85104E+01 

0 . 10000E+01 

0  *  37500E+00 
-0.5A4A7E-01 

0 ♦ 99820E+00 

"0 ♦ 4803AE-01 
-0 ♦ 99820E+00 

0 « 94837E+00 

0* AA351E-02 
-0. 10547E+01 
0*99820 E+00 

0  *  A6351E-02 
-0 ♦ 49832E-01 
-0  *  851 04E+01 

0 . 10000E+01 

0*43750E+00 
-0 ♦ 56467E-01 

0 « 99820E+00 

-0*511 09E-01 
-0 ♦ 99820E+00 

0 ♦ 94837E+00 

0.35618E-02 
-0 ♦ 10547E+01 

0 « 99820E+00 

0 • 6A351E-02 
-0  *  49832E-0 1 
-0 ♦ 85 1 04E+0 1 

0 . 10000E+01 

0  *  50000E  f 00 
-0 ♦ 24832E-02 
0.10004E+01 

-0  *  258A4E-02 
-0 • 10004E+01 

0 ♦ 99820E+00 

0 ♦ 29179E-03 
-0  * 10029E+01 

0 • 10004E+01 

0 ♦ 291 79E-03 
-0 ♦ 21914E-02 
-0 ♦ 85104E+0 1 

0 . 10000E+01 

0 . 56250E+00 
-0 ♦ 24832E-0? 

0* 10004E+01 

-0 • 275 1 8E-02 
-0.10004E+01 

0  •  99820E+00 

0  * 12A31E-03 
-0 ♦ 10029E+01 

0* 10004E+01 

0 ♦ 291 79E-03 
-0 , 2191 4E-02 
~0  *  851 04E+0 1 

0. 10000E+01 

Figure  5.  Listing  of  output  results  from  DAC  program  '<1 
execution,  case  #1. 
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0. A2500E+00 

0, 10000E+0 1 

0*35?44E -03 
“0.  lOOOOE-fOl 

0. 10004E+0I 

-0  *  49753E-04 
-0 ♦ 99960E+00 

0. IOOOOEfOl 

”0 ♦ 49753E-04 

0 ♦ 373AAE-03 
- 0  *  85104Ff 01 

0* IOOOOEfOl 

0 « 687fi0E  +  00 

0. 42342E-03 

0*  J.0000E+01 

0 «  37503E-03 
-  0. IOOOOEfO  l 

0. 10004E+01 

-0  *  27202E-04 
-0  ♦  999A0E4- 00 

0. 10000E+01 

*0  *  49753E-04 
0.37366E-03 
-0 ♦ 851 04Ef 01 

0*1 OOOOEf 01 

0.75000E+00 

0.27301E-04 

0. 10000E+0 1 

0*2A91 LF-04 
-0* 1 OOOOE+Oi 

0* IOOOOEfOl 

-0 ♦ 32079E-05 
-0  *  9?997Ef 00 

0* IOOOOEfOl 

“0  *  32079E-05 

0 « 24093E-04 
-0*85104Ef 01 

0 ♦ IOOOOEfOl 

0 ♦ 8 1250E  f 00 

0* 273C1E-04 

0. IOOOOEfOl 

0 • 28670E  04 
-0* IOOOOEfOl 

0  ♦  IOOOOEfOl 

-0* 14841E-05 
-0 . 99997E+ 00 

0*  lOOOOE-fOl 

“0 ♦ 32079E-05 

0  * 24093E-04 
-0 ♦ 85104Ef 01 

0 ♦ IOOOOEfOl 

0 « 87500E  +  00 
-0  ♦  29841E-05 

0. I0000F+05 

-0 ♦ 23842E-05 
-0* 10000E+01 

0* IOOOOEfOl 

0.35064E-06 
-0 ♦ 10000E+01 

o*iooooe+oi 

0 « 35064E-0A 
•~0  ♦  26335E-05 
“0 ♦ 85104Ef 01 

0* IOOOOEfOl 

0 ♦ 93750E+00 
0*2984  IE -05 

0* IOOOOEfOl 

~0*25034E05 
"0.1Q000E+ 01 

0* IOOOOEfOL 

0* 19791E-06 
-0*  lOOOOE-fOl 

0* IOOOOEfOl 

0 ♦ 35064E-0A 
-0.24335E-05 
-0 « 85 104Ef 01 

0  * IOOOOEfOl 

0  * 10000E+01 

0 ♦ 30515E-06 

0. IOOOOEfOl 

-0 ♦ 35763E-06 
-0* 10000F+01 

0* 10000E+01 

0 ♦ 3585AE-07 
"0 ♦ 1 OOOOEf 01 

0,  IOOOOEfOl 

0  *  35856E-07 
“0 ♦ 26929E-06 
-0 « S5104Ef 01 

0. IOOOOEfOl 

CASE  parameters: 

INTEGRATION  SCHEME  I  RUNGA-KUTTA  4TH  ORDER 
INTEGRATION  STEP  SIZE:  DT  =  0.15A25E-01 

SAMPLE  INTERVAL:  ST  -  0.12500E+00 

disturbance:  wt  -  o.iooooe+oi 

EQUATION  FOR  UNT:  UNT  -  UP  +  UD 

STEADY  STATE  OUTPUT:  X(T)  »  0.30268E-07 


Figure  5.  Listing  (concluded). 


In  addition  to  the  previous  results,  other  executions  were  performed. 

These  executions,  along  with  the  previous  resulting  execution,  may  be  summa¬ 
rized  as  in  Table  6,  and  the  resulting  graphical  plots  output  are  as  shown  in 
the  Associated  Figure} ,  column  in  the  table.  Other  executions  have  been  per¬ 
formed  with  similar  results;  however,  the  resulting  plots  appear  to  be  of  the 
same  general  shape  and  form  as  those  presented  within  this  report. 

As  can  be  seen  from  the  comparison  of  the  resultant  output  listing  and 
plots  of  Case  1  with  Case  46,  the  results  have  the  same  magnitudes  of  data 
but  are  opposite  in  sign.  Notice  that  Case  46  has  the  identical  input  data 
as  Case  1  except  for  the  opposite  signs;  they  are  mirror  images  of  each  other. 
Similar  results  have  been  obtained  for  other  cases,  such  as  Case  31  when 
compared  to  Case  32,  the  result  of  which  is  contained  in  this  report  also. 

IV.  CONCLUSIONS 

The  study  described  in  this  document  consititutes  the  first  application 
and  published  results  of  the  discrete-time  DAC  theory  to  control  systems. 

This  section  of  the  report  presents  the  conclusions  of  this  study  and  offers 
recommendations  for  further  investigations.  The  present  study  has  shown  the 
cancellation  (absorption)  of  disturbance  effects  for  the  discrete-time  DAC 
applied  to  a  control  system  for  stabilization. 

The  digital  computer  implementation  for  the  design  of  the  DAC’s  arrived  at 
for  this  example  problem  (using  the  methods  developed  in  the  appendix) 
performed,  in  each  case  of  disturbance,  those  functions  which  it  was  designed 
to  perform.  The  effects  of  the  disturbance  inputs  were  cancelled  out  by  that 
portion  of  the  controller  which  was  designed  specifically  to  handle  a  given 
waveform  mode  disturbance.  When  the  impulse  train  was  not  of  too  high  a 
frequency,  the  errors  engendered  by  the  disturbance  were  settled  out  very  well. 

A  unique  digital  computer  analysis  tool  (DDACPl — Discrete -Time  Disturbance 
Accommodating  Control  Program  1)  has  been  developed  for  implement ing  the  DAC 
control  laws,  the  equations  of  the  plant  being  controlled,  and  disturbance 
models.  Also,  a  graphical  plot  program  has  been  developed,  whereby  graphical 
plots  of  any  dependent  variable  versus  time  may  be  obtained.  Both  of  these 
programs  are  highly  ineractive  with  the  computer  user.  Additionally,  the 
graphical  plot  program  may  be  of  benefit  in  obtaining  plots  of  output  data 
from  other  programs  as  well  as  the  DAC  design  program(s) . 
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TABLE  6.  SUMMARY  OF  DAC  PROGRAM  IMPLEMENTATION  EXECUTION  RESULTS 
OBTAINED  IN  THIS  REPORT 


ASSOCIATED  FIGURES 


CASE 

NO. 

XT  AT 

T  =  0 

CWT 

AWT 

OUTPUT 

LISTING 

GRAPHICAL  PLOTS 

1 

1.0 

1.0 

0.0 

5 

6  Through  17 

2 

1.0 

1.0 

1.0 

18 

19  Through  21 

3 

1.0 

1.0 

3.0 

22 

23  Through  25 

4 

1.0 

1.0 

10.0 

26 

27  Through  29 

io1 

1.0 

1.0 

0.0 

30 

31  Through  33 

202 

1.0 

- 

- 

34 

35  Through  45 

29 

0.0 

1.0 

0.0 

46 

47  Through  48 

31 3 

1.0 

1.0 

0.0 

49 

50  Through  59 

323 

-1.0 

-1.0 

0.0 

60 

61  Through  64 

4  6 

-1.0 

-1.0 

0.0 

65 

66  Through  69 

WT  includes  a  random  noise  between  +1  input  with  a  random  number  generator 
subrout ine . 

2 

The  constant  piecewise  disturbance  was  programmed  within  the  main  discrete¬ 
time  DAC  program  as  a  function  of  time  by  the  usage  of  the  FORTRAN  IF 
statement . 

^The  program  contains  the  statement:  IF  (T  .  LT.  ST)  UNT  =  UP. 

This  statement  prevents  the  "overshoot "  in  the  plant  (system)  output 
variable  X,  and  hence,  Y. 


It  is  suggested  that  future  study  and  investigation  be  directed  to  the 
following  areas: 

1.  The  non-zero  set-point  regulator  control  problem. 

2.  The  servo-tracking  control  problem. 

3.  The  design  of  a  discrete-time  DAC  for  a  general  second-order  plant 
(system)  with  a  first-order  disturbance. 

4.  The  applications  of  discrete-time  DAC  to  a  discrete  control  problem. 
(This  would  be  beneficial  in  view  of  the  trend  toward  using  sampled-data  and 
microprocessor  techniques  in  future  designs.)  These  may  include  pointing  and 
tracking  of  designators,  gun  pointing,  autopilot  disturbance  compensation 
and  guidance  algorithm  design. 
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DAC  PROGRAM  iif  C, ASF  *2 
1  NPU I  XT  “  1.0 

FOR  EXPONENTIAL  DISTURBANCE (5) i 
INPUT  CUT  =  1*0 
INPUT  AWT  1*0 


DAC  PROGRAM  EXAMPLE  NUMBER  I  * 
OUTPUT  format: 


TIME 

XDT 

XT  ="  YT 

YNT 

UP 

UD 

UNT 

TMP:l 

XINPT 

XI  NT 

7HNT 

K 

UT 

0  *  OOOOOE+OO 

-0* 15021E+02 

0* lOOOOE+Ol 

0* 10000E+01 

0 • 85104E+01 

-0 ♦ 85104E+01 

-0* 17021E+0? 

-0 ♦ 80000E+01 

0.90A52E+01 

0* OOOOOE+OO 

0 . SOOOOE+O 1 

~0 ♦ 85104E+01 

0* 10000E+01 

0, t 0000E+01 

-0. 15099E+00 

0.37253E-07 

0*37253E-07 

•0.31704E-06 

-0 ♦ 28693E+0I 

-0.28 693E+01 

-0 ♦ 29B02E-06 

0 ♦ 30563E+01 

0 ♦ 2A972F+01 

0  *  2A972E+01 

-0 ♦ 85104E+01 

0 ♦ 27183E+01 

CASF  PARAMETERS : 

INTEGRATION  SCHEME:  RUNRA-KTJTTA  4TH  ORDER 
INTEGRATION  STEP  SIZE:  DT  **  0  *  :l  5A25E-01 

SAMPLE  INTERVAL:  ST  ~  0  *  1 2500E+00 

DISTURBANCE:  UT  -  0.27183E+01 

EQUATION  EOR  UNT :  UMT  UP  +  UD 
STEADY  STATE  OUTPUT:  X<T>  -  -0.23591E-02 


Figure  18.  Output  listing  (condensed)  of  results,  Case  #2. 


38 


Figure  19.  Plot  No.  1,  DAC  progran  if  1,  Case  ftl. 
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DAC  PROGRAM  #1*  CASE  t 3 
INPUT  XT  ~  1 ♦ 0 

FOR  EXPONENTIAL  DISTURBANCE ( S ) : 
INPUT  CWT  -  1*0 
INPUT  AUT  *  3*0 


DAC  PROGRAM  EXAMPLE  NUMBER  1* 
OUTPUT  format: 


TIME 

UP 

XINPT 

3  ♦  OOOOOE-f  00 
D ♦ 85  L04E+01 
D ♦ 13155E  f 02 

!K10000E+01 
!>*19022E-0<S 
3 ♦ 28547E+ 02 


XDT 

UD 

XINT 

-0* 17438E+02 
-0  * 10928E+02 
0 ♦ OOOOOE+OO 

-0*36272E+01 
“0 • 23713E+02 
0 ♦ 19620E+ 02 


XT  a  YT 

UNT 

ZHNT 

0* 10000E+01 
-0*  L9438E+02 
Q.90416E+01 

-0  *  22352E-07 
-0  *  23713E+02 
0* 19620Ef02 


YNT 

TMP:L 

K 

0 ♦ 10000E+01 
-0 ♦ 90416E+ 01 
-0*  85104E+01 

-0 ♦ 22352E-07 
0 « 20210E-06 
-0*85104E+01 


0 ♦ 10000E+01 


0.200Q6E+0 2 


CASE  parameters: 

INTEGRATION  SCHEME:  RUNGA-KUTTA  4TH  ORDER 
INTEGRATION  STEP  SIZE:  DT  »  0  ♦  iri625E~01 

SAMPLE  INTERVAL:  ST  «  0,12500E+00 

DISTURBANCE:  WT  *  0*20086E+02 

EQUATION  FOR  UNT:  UNT  a  UP  4  UD 
STEADY  STATE  OUTPUT:  X(T)  -  -0 ♦ 56675E-01 


Figure  22.  Outnut  listing  (condensed)  of  results.  Case  #3. 
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r» AC  PROGRAM  *1  It  GASP  t4 
INPUT  XT  -  1.0 

FOR  EXPONENT I Al  DISTURBANCE < S ) : 
INPUT  CWT  »  1.0 
INPUT  AWT  -  10.0 


DAC  PROGRAM  EXAMPLE  NUMRFR  I 
OUTPUT  format: 


TTMF 

UP 

XINPT 

0 ♦ OOOOOF  f 00 
”0 ♦ 851. 04E+01 
0.4A514E+02 

0. lOOOOE+Ol 
-0.324A5E-03 
0*7105  IE +05 


-0.32724E+0 2 
-0.2A214E+02 
0 ♦ OOOOOEf 00 

-0*1 801 6F+05 
-0.40043E+05 
0 • 20357E+05 


0. LOOOOFfOl 
-0 . 34724E4-02 
0, 13326E+02 

0.38147E-04 
-0 . 40043E+05 
0  *  20357E+05 


0 ♦ lOOOOE+Ol 
-0 ♦ 1332AE+02 
-0 ♦ 85104E+01 

0 ♦ 381 47E-04 
-0.5083AE-03 
-0 . 85104E+01 


0. 10000F+01 


0  *  22026F+05 


CASE  parameters: 

INTEGRATION  SCHEME:  RUNGA-KUTTA  4TH  ORDER 
INTEGRATION  STEP  SIZE:  »T  *  0.15A25E-01 

SAMPLE  INTERVAL:  ST  -  0*12500Ff00 

disturbance:  wt  *  0.22026E+05 

EQUATION  FOR  UNTi  UNT  *  UP  +  UD 

STEADY  STATE  OUTPUT:  X(T>  *  -0.28150E+03 


Figure  26.  Output  listing  (condensed)  of  results  for  Case  #4. 
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DAC  PROGRAM  *1,  CASE  *10 
INPUT  XT  •  1.0 

FOR  EXPONENTIAL  DISTURBANCES )» 
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dH«4 

d  h 
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UJ  UJ  Ul 

nr-  ▼ 
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mood 

UJOOd 
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oo<* 

in  ^  in 

in  ^  in 

00  Ul  00 

00  Ul  00 

ddd 

ddd 
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Figure  30.  Output  listing  of  results  for  Case  *‘10  (continued'). 


0.93750E+00  0.10696E+01  0.91363E-02  -0.25347E-01 

0.21571E+00  -0.75639E+00  -9.54968E+90  0.19936E+00  0.60116E+00 

0.75639E+00  0.94676E+00  0.75639E+00  -0.85104E+01  0.16012E+01 


m  ai 
♦  i 

yui 
(0  AI 
r-  h 
oo  t* 
OMU 


i  ♦  ♦ 
UllAlUl 
tt(U^ 

am# 

aau) 

maa 


i  i 


oowr- 
a^a 
~«coin 
am# 
m  'H 


•  i 


H  ®  H 
### 


®  H  Is- 

#U)# 

#mu> 
am# 
«^m^ 
♦  •  • 
### 
i 


g< 


10 

OJ 

a  < 

u>  < 


C  *4)  O  • 
h##(ud# 
fc  a# 

3  U>  I  f 
w  •  oiiai  ■ 
1  ra  a. 


H  H# 

tty  • 

ft)  I  w 

H  x 
•  1- 

•  00  O 

M 

♦  ♦  ♦ 

«•  ML 

3  ft- 

UIUIUI 

uj  in  in 

oo  r-  oi 

■  a. 

m#m 

••  y  a.  «* 

i- 1- 

▼  inu> 

in  z  uj  fj 

>-Z3 

m  #  in 

ttOI-tt 

3  3  0 

•  •  • 

uj  in  in  3 

0-  tt 

M  ttUl 

#  ®  ® 

UZZUIUOK 

OOI-OIl# 

nnr  r  k* 
•-hnara 

<E  <L  0#  O 
Q.aayocM> 
0(9  J3hQ 
U1UJUICL  I-  «  a 

55tt5w3w 

«CZ<EHOh 

OHM#oy  <n 


CL 

O 

I- 

in 


CL 

<n 

o 


52 


Fipure  30.  Output  listing  of  results  for  Case  #10  (concluded). 


UN  SDACP? 


BAG  PROGRAM  fir  CASF  t?.0 
INPUT  XT  -  1.0 


nor  PROGRAM  EXAMPLE  pHIMBFR  I  . 
OUTPUT  format: 


TTMF 

UP 

XINPT 

0  • OOOOOEf 00 
-0.85104E+01 
0  ♦  7f5:l  04E+01 

0.42500E-01 
-0.95104E+01 
0 . 75 1 04E+01 

0.12500E+00 
0.724A7F.+01 
0.  U153E  +  01 

0 . 18750E+00 
0.724<67E-f0t 
0. 11153E+01 

0  *  25000E+00 
0, 18912E+00 
0 ♦ 94837E+00 

0 • 31250E+00 
0, 18912E+00 
0 ♦ 94837E+00 

0.37500E+00 
0 . 32105E+01 
-0.18849E+01 

0 ♦ 43750E+00 
0 ♦ 32105E+01 
-0 ♦ 18849E+01 

0 ♦ 50000E+00 
0.35272E+01 
-0.49977E+01 

0.56250E+00 
0 . 35272E+ 0 1 
-  0 . 49977E+0 1 

0 . 62500E+ 00 
0  *  31607E-0 J 
-0.5025AF-f0:t 


XriT 

(JO 

XT  NT 

-  0 ♦ 1 402 IE  FO? 
-0.75I04E+01 
0 . OOOOOEf 00 

-0. 1 491 8E 4*02 
-0.75104E+01 
0. OOOOOEf 00 

0 . A7799E+01 
-0. 11153E+01 
0.75104E+01 

0.6A817F+0J 
-0  ♦  .1 1 1 53E+01 
0 ♦ 75 104E+01 

0 ♦ 21853E+00 
~0.94837E-fOO 
0.1 1 153E+01 

-0 ♦ 57A75E+ 01 
-0.94837E+00 
0.1 1 153E+01 

-0 ♦ 28182E+00 
0 ♦ 18849E+01 
0.94837E+00 

-0 ♦ 29986E+00 
0 ♦ 1 8849E+01 
0  ♦ 94837E+00 

0*31 I05E+01 
0  * 49977ET 01 
-0  ♦  18849E+01 

0 ♦ 33095E+01 
0.49977ET01 
-0  ♦  18849E+01 

0.53461E—01 
0 .  5025AF+01 
-0 ♦ 49977E+C1 


XT  ™  YT 

un  r 

ZHNT 

0. 10000F+01 
0, I A021ET02 
0.751 04E+0J 

0. J 0294ET 00 
-0. IA021E+02 
0.75104E+01 

-0 ♦ 851 50E  fOO 
0.61314E-F01 
0.11 153E+01 

-0  ♦  44972ET 00 
0.A1314E+01 
0.11 .153E+01 

-0 ♦ 22222E -01 
“0 ♦ 75925EL 00 
0 . 94837E+ 00 

-0  ♦  824Q8E  ~02 
-0 ♦ 75925E+00 
0 ♦ 94837E+00 

-0.37725E+00 
0.50954E+01 
-0. 18849E+01 

-0 ♦ 39528E+ 00 
0.50951E+01 
-0. 18849E+01 

-0. 4  L446E+00 
0.85249E+01 
-0  *  49977E4- 01 

-0  *  21545E+00 
O.85249E+01 
-0 ♦ 49977E+01 

-0 . 371 33E-0? 
0 . 50572E+01 
-0  -  FI025AE+01 


YNT 

TMPi 

I 

0. JOOOOE+Ol 
-0 ♦ 75104E+01 
-0.85 104E+01 

0. 10000E+01 
•0 , 75104E+01 
-0.85104E+01 

-0.85150E+00 
0 . 63951E+01 
-0. 851O4E+01 

“0 . 85150E+00 
0.A3951E+01 
— 0.85104E+01 

-  0 , 22222E-01 
0.16690E+00 
-0. 85104E+01 

”0  ♦  22222E-01. 

0 ♦ 1A690E+00 
-0  *  85104ET01 

~0. 37725E+00 
0.28333E+01 
-0.85104E+01 

-0.37725E+00 
0. 28333E+01 
0.85104ET01 

—0 • 4 1 446E+00 
0 . 31 128E+01 
-0.85104E+01 

”0.414  46E+00 
0.3M28E+01 
”0 ♦ 85104E+01 

-0 . 37133E-02 
0.27889E-01 
-0.85104E+01 


Figure  34.  Output  listing  of  results  for  Case  #20. 
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f  * 


0. 10000E+01 


0.10000E+0 1 


0.1 OOOOE-f  01 


0. 10000E+01 


0. tOOOOE+Ot 


”0 . 50000E+ 01 


-0 « 50000E+01 


”0 • 50000E+ 01 


-0.50000E+01 


”0 . 50000E+0 1 


-0.50000F+0 I 


0.AB7S0F+00 

0.31<S0:?E-0i 

-0.50756F+0I 

0  *  5A981E—01 

0  * 5025AE+01 
-0 . 49977F+01 

~0  * 292B9F -03 
0. 5O“i72r+01 
-  0  ♦  502f»AF+01 

-0 ♦ 37133E-02 

0 « 27889F-01 
-0  *  85104E+01 

-0.50000E+01 

0 . 75000E+00 
-0.20fv'14F  +  01 
-0.37i?5F+01 

0* 13995F+01 
0*321 25E  f 01 
-0  ♦  5025AE4- 01 

0*241 4 OF  f 00 
0* 1 1 58 lE+Ot 
~0- 32175F+01 

0*241 40E+00 
-0* 19130E+01 
-0 . 85104F+01 

0 • 00000F+00 

0  *  91 250E+00 
-0.20544E+01 
-0*32125F+01 

0* 14890F+01 

0 , 32125F  f 0 1 
-0 ♦ 5025AF+01 

0 . 33094E+00 
0* 1 J581E+01 
-0 . 321 25F+01 

0*24 1 40E+00 
”0 ♦ 1 81 30E+01 
-0*85104F+01 

0 ♦ 00000F+00 

0  ♦  87500E+00 
-0 ♦ 36272E+01 
-0, 1150AE-01 

-0 • 31 895E+ 01 
0*1J  50AF-01 
-0 ♦ 32125F+01 

0. 42A21E+00 
-0  *  3A157E+01 
-0*1 150AE-01 

0. 42421E+00 
-0 ♦ 32010E+01 
-0  *  851 04E+01 

0 . OOOOOE+OO 

0.93750E+00 
-0 ♦ 3A272E+ 01 
-0. 11506E-01 

-0 . 3393AF+01 

0. 11506F-01 
-0.32125E+01 

0 , 22215F+00 
- 0.3A157E+01 
-0* 1150AE-01 

0 ♦ 42A21 E+00 
“0 • 32010F+01 
-0 « 851 04E+01 

0 ♦ OOOOOE+OO 

0* 10000F+01 
-0 . 42731E-01 

0  *  26222E-01 

0 ♦ 993A1E+ 01 
-0 ♦ 2A222E-01 
-0* li50AE~01 

0 ♦ 50234F-02 
-0 . A8973E-01 

0 . 26222F-01 

0 . 50234F-02 
“0  ♦  3772BF--01 

♦ 85104E+0 1 

0* 10000E+02 

0. 10A25E+01 
-0*42751E~01 

0 ♦ 2 A222E-01 

0* 10572E+02 
-0 . 2A222E -01 
-0.11306E-01 

0  *  64073F  f 00 
“0  *  A8973E-01 

0  *  26222F-0 1 

0 ♦ 50234F-02 
-0  *  37728E-01 
-0  *  85104E+01 

0 ♦ 10000E+02 

0  ♦  1 1250E+01 
-0.98795E+01 

0  *  87448E+01 

-0* 174A3E+02 
-0 « 87448E+ 01 

0 ♦ 2A222E-01 

0* 11A09E+01 
-0 ♦ 18A24E+02 

0 ♦ 87448E+ 01 

0.11 A09F+0 1 
-0  *  87 1  SAF'f  01 
-0  *  85104F+01 

O.OOOOOE+OO 

0*1 1875E+01 
-0.98795E+01 

0  *  87448E+01 

-0*1 858 1E+02 
~0 • S7448E+01 

0 ♦ 26222E-01 

0  *  43555F-01 
-0*  18A24E*f02 

0  * 87448F+01 

0*  11A09F+01 
-0 . B718AE+01 
”0  *  85104E+01 

0  4  OOOOOE+OO 

0* 12500E+01 
0*97465E+01 

0* 14357E+00 

0 ♦ 84577E+01 
-0* 14357E+00 

0 ♦ 87448F+01 

-0*  U452E+01 

0 ♦ 9A029E+01 

0  *  1 4357E+ 00 

-0* 11452E+01 

0 « 8601 2E+01 
-0  *  35104E+01 

0 . OOOOOE+OO 

0* 13125E+01 
0*974A5E+0t 

0. 14357E+00 

0  *89988E+01 
-0 « 14357E+00 

0*  87448E+01 

-0 ♦ A0412E+00 
0*94029E+01 

0. 14357F+00 

-0 ♦ 1 1 452E+01 

0*  86012E+01 
-0 ♦ 85104E+01 

O.OOOOOE+OO 

0 ♦ 13750E+01 
0*241 48E+00 
-0 ♦ &9533F-01 

0 • 28264E+ 00 

0 ♦ A9533E-01 

0* 1 4357E+00 

-0 ♦ 28374F-01 

0  ♦  31  lOlF-fOO 
-0*69533F~01 

-0  *  28374E-01 

0 ♦ 21310E+00 
-0 ♦ 85104F+01 

0. OOOOOE+OO 

0 « 14375E+01 

0  *  24 1 48F+00 
-0.69533E-01 

0 « 30072F+00 

0 ♦ 69333E-01 
0*1 4357F+00 

-0  ♦ :!  0291E-01 

0 ♦ 31 101F+00 
-0 ♦ A9533F-01 

-0  *  28374E-01 
0.21310E+00 
-0 ♦ 85104E+01 

O.OOOOOE+OO 

0, 15000E+01 
-0.7A158F-01 
-0 . 2323AF -02 

-0.A488SE-01 

0 . 2323AF -0? 
-0.  A9533E-01 

0 . 894R8F-02 
”0  *  73S34F  ~01 
-0.2323AE-02 

0  *  89488E -02 
~0*A7209F-0t 
-0 ♦ 85104F+01 

0 . OOOOOF+OO 

F^urt*  34.  Output 

listing  oT  results 

for  Case  #20  (continued). 
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sdacp: 


CASH  parameters: 

INTEGRATION  SCHEME:  RIJNGA-KUTTA  4TH  ORDER 
INTEGRATION  STEP  SIZE:  DT  =-  0 ♦ 15A25E-01 

SAMPLE  INTERVALS  ST  ^  <K12500Ei00 

disturbance:  wt  -  o*oooooek>o 

EQUATION  FOR  I.JNT:  UNT  ~  UP  +  IJD 

STEADY  STATE  OUTPUT:  X<T)  -  0*79349E-02 


—  STOP 


;! 

'  j 
-  i 


i 


Figure  34.  Output  Listing  of  results  for  Case  #20  (concluded). 

> 
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PLOT  NO.  7;  DAC  PROGRAM  #1,  CASE  #8G 


Figure  41.  Plot  No.  7,  DAC  program  #1,  Case  #20. 


Figure  44.  Plot  No.  10,  RAC  program  #1,  Case  #20. 
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I* AC  PROGRAM  #1, 
INPUT  XT  -  0.0 
FOR  EXPONENTIAL 
INPUT  CUT  -  1*0 
INPUT  AUT  -=  0.0 


CASE  *29 
DISTURBANCE'S) ! 


BAC  PROGRAM  EXAMPLE  NUMBER  I . 

output  format: 


YNT 

TMPl 

K 

0  #  OOOOOE+OO 
0 . OOOOOE+OO 
“0 ♦ 85104E+01 

~0. 13970E-07 
0.10492E-06 
-0 ♦ 85 104E+01 


CASE  parameters: 

INTEGRATION  SCHEME :  RUNGA-KUTTA  -1TH  ORDER 
INTEGRATION  STEP  size:  DT  -  0.15625F-0I 

SAMPLE  INTERVAL:  ST  -  0.12500E+00 

DISTURBANCE :  WT  »  0.10000E+01 

EOUAT  ION  FOR  UNT :  UNT  ~  UP  +  IJD 
STEADY  STATE  OUTPUT:  X(T)  a  -0.12107E-07 


TTMF 

UP 

XINPT 

0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 

0. 10000E+01 
0*1 1BR9F-06 
C. lOOOOE+Ot 


XDT 

UD 

XINT 

0. ) OOOOF+O I 
0. OOOOOE+OO 
0. OOOOOE+OO 

0 . i 1921E-0A 
-0. 10000E+01 
O. lOOOOE+Ol 


XT  *  YT 

UNT 

7HNT 

0  #  OOOOOE  +  OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 

-0. 13970E-07 
-0.10000E+01 
0. 10000E+01 


TTO  —  STOP 


Figure  46.  Output  listing  (condensed)  of  results  for  Case  #29. 


UT 


0. 10000E+OI 


0.1 OOOOE+Ol 
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DAC  PROGRAM  t  U  CASE  *31 
INPUT  XT  =s  1,0 

FOR  EXPONENTIAL  DISTURBANCE C S ) J 
INPUT  CUT  »  1.0 
INPUT  AUT  «  0.0 


DAC  PROGRAM  EXAMPLE  NUMBER  1. 
OUTPUT  format: 


TIME 

XDT 

XT  «  YT 

YNT 

UP 

UD 

UNT 

TMP1 

XINPT 

XINT 

7HNT 

K 

WT 

0 ♦ OOOOOE+OO 

-0.A5104E+01 

0. 10000E+01 

0 ♦ 10000E+01 

-0 ♦ 85104E+01 

-0 ♦ 75104E+01 

“0 ♦ 851 04E+01 

-0 ♦ 75104E+01 

0. OOOOOE+OO 

0. OOOOOE+OO 

0 . 75104E+01 

-0 ♦ 85 104E+01 

0 ♦ 10000E+01 

0. 10000E+01 

0. OOOOOE+OO 

0.93132E-09 

0 ♦ 931 32E-09 

-0.79259E-08 

-0. 10000E+01 

“0. 10000E+01 

-0.69946E-08 

0 • 10000E+01 

0 * 10000E+01 

0 ♦ 10000E+01 

-0.85104E+01 

0  * 10000E+01 

case  parameters: 

INTEGRATION  SCHEME:  RUNGA-KUTTA  4TH  ORDER 
INTEGRATION  STEP  SIZE:  DT  *  0.15625E-02 

SAMPLE  INTERVAL:  ST  «  0.12500E+00 

disturbance:  wt  ®  o.iooooe+oi 

EQUATION  FOR  UNTJ  UNT  »  UP  +  UD 

STEADY  STATE  OUTPUT:  X(T)  «■  (K93132E-09 


Figure  49.  Output  listing  (condensed)  of  results  for  Case  *31. 


! 
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PLOT  NO.  7)  DAC  PROGRAM  *1,  CASE  t31 


Figure  36.  Plot  No.  7,  DAC  program  // 1,  Case  1*31. 
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PLOT  NO.  9;  DAC  PROGRAM  #1,  CASE  «31 


PLOT  NO.  1}  DAC  PROGRAM  *1,  CASE  #32 


#V 


Figure  61 .  Plot  No.  1,  DAC  program  ^1,  Case  If 32 
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T'AC  PROGRAM  *1»  CASE  #46 
INPUT  XT  ~  -1.0 

FOR  EXPONENTIAL  BISTURBANCE(S)  J 
INPUT  CUT  ~  -1.0 
INPUT  AUT  -  0.0 


DAT.  PROGRAM  EXAMPLE  NUMBER  1  ♦ 
OUTPUT  format: 


1 1  nr 

UP 

XINPT 

XBT 

UD 

XINT 

XT  -  YT 

UNT 

2HNT 

YNT 

TMP1 

K 

UT 

0 ♦ 00 000 FI 00 
0.R3 L04F10 1 
-0 . 7fi  1  0  4F+0  1 

0.1 4021E+02 

0 . 75 104E+01 

0 . OOOOOE+OO 

-0. JOOOOE+Ol 
0.16021E+02 
-0 . 751 04E+01 

-0. 10000E+01 
0.75104F+01 
-0 . 85104E+01 

-0 . 10000E+01 

0 . A7500r-0 » 
0.851 04F+0J 
-  0.751 04Fio:l. 

0. 14918E+02 

0 . 751 04E+01 

0. OOOOOF+OO 

-0. 10294F+00 

0. 1A021E+02 
-0 . 751 04F+01 

-0 . 10000E+01 
0.75104E+01 
-0 . 851 04E+01 

-0. 10000E+01 

o.  i  oo 

-o.7?4<v:7r  foi 
-o.  t  i  i?»3i:+o x 

-0.62799E+01 

0 , 1 1 153E+0 1 
-0.75104E+C1 

0.85150F+00 
-0.6131 4E+0.1 
-0. 1 1 153E+01 

0 . 85150E+00 
-0 . 63951E+01 
-0 . 851 04E+01 

-0. 10000E+01 

0. I8750E+00 

0. 7?4*7E+0J 
-  0.1  U53E+01 

-0.6681 7F+0I 
0.11 153E+0J 
-0 . 75104E+01 

0 . 44972E +00 
-0.6131 4E+01 
-0. 11153E+01 

0.85150E+00 
-0.63951E+01 
-0. B5104F+01 

-0 . 10000E+01 

0.2 50 OOF  100 
-0..I8912E100 
-0.94837E+00 

0.21853E+00 

0 . 94837E+00 
-0. 11153E+01 

0 . 22222F-01 

0 . 75925E+00 
-0 . 94S37E+00 

0.22222E-01 
-0 . 16690E+00 
-0 . 85104E+01 

-0 . 10000E+01 

0.31250F+00 
-0. 1891 2E +00 
-0.94837E100 

-0.23251E+00 

0 . 94837E+00 
-0.11 153E+01 

0 . 82408E-02 

0 . 75925E+00 
-0 . 94837E+00 

0.22222E-01 
-0. I6690E+00 
~0 . 85104E+01 

-0 . 10000E+01 

0 ♦ 37500E+00 
0.5A4A7F-01 
-0.99820F100 

0 . 48036E-01 

0 , 99820E+00 
-0 . 94837E+00 

-0.66 351E-02 

0. 10547E+01 
-0 . 99820E+00 

-0.66351E-02 

0, 49832E-01 
-0.85104E+01 

-0 . 10000E+01 

0. 43750F100 
0.5A4A7F-01 
-0 . 99820F100 

0.5  109F-01 
0.6 ;820F+00 
-0 . 94837E+00 

-0.35618E-02 

0. 10547E+01 
-0 . 99820E+00 

-0 . 66351F-02 
0.49832E-01 
-0 . 85104E+01 

-0. 10000E+01 

0 ♦ 50000F4 00 
0.24832E-02 
-0. 10004E+01 

0.25864E-02 

0. I0004E+01 
-0. 99820E+00 

-0.29179E-03 
0.10029E+01 
-0. 10004E+01 

-0.29179E-03 

0 . 21 91 4E-02 
-0 . 85104E+01 

-0.10000E+01 

O.5A250F+00 
0.24832E--02 
-0. 10004E+0L 

0 . 27518E-02 

0. 10004E+01 
-0 . 99820E+00 

-0 . 12631E-03 

0* 10029E+01 
-0. 10004E+01 

-0 . 29179E-03 

0 . 21914E-02 
-0 . 85104E+01 

-0. 10000E+01 

Figure  65.  Output  listing  of  results  for  Case  #46. 
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